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Abstract 

We study a two-component asymmetric simple exclusion process (ASEP) 
that is equivalent to the ASEP with second-class particles. We prove self¬ 
duality with respect to a family of duality functions which are shown to arise 
from the reversible measures of the process and the symmetry of the generator 
under the quantum algebra Uq[Ql^]. We construct all invariant measures in 
explicit form and discuss some of their properties. We also prove a sum rule 
for the duality functions. 
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1 Introduction 


We consider an asymmetric simple exclusion process with two species of particles 
on the one-dimensional hnite lattice A = {—L -|- 1,..., L}. Its Markovian dynamics 
can be described informally as follows. Each site i can be either empty (denoted by 
0) or occupied by at most one particle of type A or of type B. Thus we have local 
occupation numbers ri{k) G {A,0,i?}. We dehne the bonds {k,k -|- 1) of A where 
—L + l<k<L — 1. Each bond carries a clock i which rings independently of all 
other clocks after an exponentially distributed random time with parameter where 
Tfc = r if {T]{k),r]{k -h 1)) G {(A, 0), (0, B), (A, B)} and Tk = Hi {ri{k),r]{k + 1)) G 
{(0, A), (i?,0), {B,A)}. When the clock rings the particle occupation variables are 
interchanged and the clock acquires the corresponding new parameter. Symbolically 
this process can be presented by the nearest neighbour particle jumps 

AO ^ OA ') 

OB —)■ BO > with rate r 
AB ^ BA J 

OA ^ AO ') 

BO —)■ OB > with rate i 
BA^AB J 

We have reflecting boundary conditions, which means that no jumps from the left 
boundary site —L-l-1 to the left and no jumps from the right boundary site L to the 
right are allowed. We shall assume partially asymmetric hopping, i.e., 0 <r,i < oo. 
By interchanging the i?-particles and vacancies this process turns into the ASEP 
with second-class particles [1] . We choose an even number of lattice sites exclusively 
for the sake of convenience of notation. 

The objective of this work is to construct for the hnite lattice in explicit form all 
reversible measures and to prove self-duality with respect to a family of duality func¬ 
tions that allows for the computation of expectations of the many-particle system 
in terms of transition probabilities of the same process with only a small number 
of particles. It will transpire that this property, analogous to the well-known self¬ 
duality of the simple symmetric exclusion process [2], arises from the fact proved in 
[3j that the generator of this process commutes with a set of matrices which form a 
representation of the quantum algebra f/q[0[(3)], which is the g-deformed universal 
enveloping algebra of the Lie algebra 0l(3) dehned below (([22]) - fl26|) j. 

The idea of deriving of duality relations from the representation matrices of 
a non-abelian symmetry algebra of a generator of a Markov process goes back to 
Schiitz and Sandow [1] where this strategy was applied to the symmetric partial 
exclusion process on arbitrary lattices. This is an interacting particle system with 
a SU{2) symmetry where each lattice site can be occupied by at most a hnite 
number of particles. Next this symmetry approach was extended to prove self¬ 
duality of the asymmetric simple exclusion process [Sj, which is symmetric under 
the action of quantum algebra f/g[ 0 l( 2 )] and which is an integrable model solvable 
by Bethe ansatz. The self-duality together with the integrability was used in [6j to 


( 1 ) 

( 2 ) 
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study the time evolution of shock measures and in to study current moments. 
By mapping the ASEP to a lattice model of interface grwoth the duality function 
can be interpreted as a lattice Cole-Hopf transformation [5] and is therefore yields 
information on KPZ interface growth and the moments of the partition function of 
a directed polymer [8]. 

The idea of using symmetries of the generator to obtain duality functions was 
employed again by Giardina et. ah [9] to study heat conduction in the KMP model 
with SU{1,1) symmetry and subsequently extended to other interacting particle 
systems, including particle systems without conservation of particle number m 
[m [El [13]. Recently the f/q[gl(2)] symmetry was extended to the non-integrable 
asymmetric generalization of the S't/(2)-symmetric partial exclusion process [IT] . 
Duality relations for new integrable models that can be solved by Bethe ansatz and 
related methods were studied very recently in [151 [16] . 

Here we prove self-duality for the two-component ASEP mentioned above whose 
symmetry algebra f/g[gl(3)] is larger than SU{2), SU{1, 1) or their g-deformations. 
We shall consider only hnite systems, the construction and characterization of the 
properties of the process on Z is out of the scope of this work. The main novel 
feature is the presence of more than one conserved species of particles. This leads 
to interesting non-local properties of the duality functions and, through the inte- 
grability of the model, to the possibility of applications in the inhnite volume limit 
employing exact computations along the lines of [T71 [T8] . 

The paper is structured as follows. In Sec. 2 we define the process and mention 
two results obtained in recent work [3] that will be used here. In Sec. 3 we state 
the main results of the present work. Sec. 4 is included for self-containedness. We 
describe some tools from linear algebra [191 [20] used in the proofs, which are con¬ 
venient, but not widely known in the probabilistic treatment of interacting particle 
systems. In Sec. 5 we present the proofs of our results. 


2 The two-component ASEP 

We define the process, introduce notation, and mention some results used in the 
proofs. 

2.1 State space and configurations 

It is convenient to introduce ternary local state variables rjlk) G S where S = 
{0,1,2}. We say that 0 represents occupation of a site a particle of type A, 1 
represents a vacant site and 2 represents occupation by a particle of type B. Thus 
a configuration is denoted by 77 = {ri{—L + 1),... ,ri{L)} G S^'^. We call this 
characterization of a conhguration the occupation variable representation. We shall 
repeatedly consider configurations with a hxed number N particles of type A and 
M particles of type B. We denote configurations with this property by 77 atm and 
the set of all such conhgurations by 
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Equivalently we can specify a configuration rj uniquely by indicating the particle 
positions z on the lattice and write 


z = {x,y} 


(3) 


with 

X ;= {xi : ri{xi) = 0}, y ;= {yi : r]{yi) = 2} 

We call this the position representation. 

Throughout this work we use the Kronecker-symbol dehned by 


f 1 ii a = (3 

\ 0 else 


(4) 

(5) 


for a, (3 from any set. We also introduce for k,l ^ A 


and the indicator function 


lA(r7) 


1 ii r] & A 
0 else 


( 6 ) 

(7) 


for subsets A C S^'^. Some other functions of the conhgurations will play a role in 
our treatment: 


Definition 2.1 For 1 < k < L we define the local permutation 

^ {r^(-L + l),...T]{k-l),r]{k + l),ri{k),T]{k + 2),...,r]{L)} =: (8) 

Definition 2.2 We define local occupation number variables 

Otkidl) ■ ^kijfi) ■ drj{k'),2 (9) 

and the global particle and vacancy numbers 

L L L 

^i'n)= Yl M{r))= Y V{r))= Y (10) 

k=—L-\-l k=—L+\ k=—L+\ 

The argument of the local occupation number variables will be suppressed through¬ 
out this paper, but not the argument of the global particle and vacancy numbers. 
We note, for z = r], the trivial but frequently used identities 


N{r]) = N{z) = |x|. 

M{r]) = M(z) = |y|. 

(11) 

V(z) 

M(z) 


i=l 

bk = 'Y^ byi,k- 
i=l 

(12) 
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Definition 2.3 For configurations z = {x, y} we define the number Nk{z) of A- 
particles to the left of site A; G A and analogously the number Mk{z) of B-particles 
to the left of site /c G A 

fc-l N{z) k-1 k-1 M(z) k-1 

jvt(z):= Yi = E E E ‘‘“E E (13) 

Z=—L+l i=l l=—L+l l=—L+l i=l l=—L+l 

and 

Ak{z) := 2Nk{z) - iV(z), Bk{z) := 2Mfc(z) - M(z). (14) 

Notice that the functions N{-), Nk{-),ak{-), Ak{-) depend only on the x-coordinates 
(positions of the A-particles) of a configuration z, while M{-), Mk{-),bk{-), Bk{-) 
depend only on the ^-coordinates. 


2.2 Definition of the two-component ASEP 

Recalling the definitions (IH]) and flT^ the two-component ASEP rj^ described infor¬ 
mally in the introduction is defined by the generator 

L-l 

£/(>?)= E (15) 

k=—L-\-l 

with the local hopping rates 

y^fcfc+i(^) ^ ^ {akVk+1 + Vkbk+i + ttkbk+i) + i {vkak+i -h bkVk+i + bkOk+i) (16) 


for a transition from a configuration 77 to a configuration ry' with transition rate 

L-l 

w{r] ^ rj') = ^ (17) 

k=—L-\-l 


It will turn out to be convenient to introduce the asymmetry parameter q and time- 
scale factor w 


q = 


w = 


y/rl. 


(18) 


The time scale will play no significant role below. 

The general form of the evolution equation of a Markov chain with state space 
and transition rates w{rf ^ r}') for a transition from a configuration 77 G to a 
configuration 77 ' G is 


^/(^) = ^ h')[/(h') - /(h)] 


(19) 


where the prime at the summation indicates the absence of the term rj' = rj. We 
define the transition matrix H of the process by the matrix elements 


H 


v'v 


—w{r] —)■ 77 ') 77 7 ^ 77 ' 


( 20 ) 
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The defining equation (fT^ then becomes 

£/(»?) = -E( 21 ) 

rj'GO 

The r.h.s. of (12T]) represents the multiplication of the matrix H with a vector whose 
components are /(r;') in the canonical basis. In slight abuse of language we shall 
call also H the generator of the process. Below we shall construct H for the two- 
component exclusion process in a judiciously chosen basis. 


2.3 The quantum algebra UMn)] 

For the Lie algebra fll(n) the quantum algebra Uq[Q{{n)] is the associative algebra 
over C generated by i = 1 , ..., n and Xf, i = 1,..., n — 1 with the relations 

[2I1E2] 


[ Lj, Lij] — 0 

^-1 . - (WiLri)-2 

1 ] = %-- 

and, for 1 < f, j < n — 1, the quadratic and cubic Serre relations 

[xyxp = o K-iljii, 

(XffXf - I2],xtxfxt + Xf(Xff = 0 \i-j 

Here the symmetric g-number is defined by 



qx _ q X 

9-9"' 


1 . 


( 22 ) 

(23) 

(24) 


(25) 

(26) 


(27) 


for 7 I 0 and x e C. (Notice the replacement 9 ^ —t 9 that we made in the 

definitions of [ 22 ]-) The notion of symmetry of the generator under the action of 
the algebra means that there exist representation matrices and L, of the algebra 
that all commute with the transition matrix H of the process, i.e., 


l.f/,K±] = [ff,L.| = 0. 


(28) 


For the present case n = 3 these representation matrices, given in fllOSp and (11171) . 
were constructed in [3]. 


3 Results 

In order to state the first main result we first define the g-factorial 

n>l 


(29) 
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and the g-multinomial coefficients 

|A'],! 


CkW = 


|]V],![A--]V], 


, CK(iV,M) = 


[K] 


<?• 




(30) 


Theorem 3.1 The two-component exclusion process m restricted to the subset 
^ particles of type A and M particles of type B has the unique invariant 

measure 

with the reversible measure 

Tlir]) = (^k^-L+d‘^^~^)Ak—bk)+'Bk=-L+lTil^-L + d°-i^k+l—ho-k+l) (32) 

and the normalization factor 

Z2l{N,M) = C2l{N,M). (33) 


We shall call these invariant measures, characterized in the following theorem, 
the canonical equilibrium distributions of the process. Particle number conservation 
yields the following corollary. 


Corollary 3.2 The convex combinations 


2L 2L-N 


QlA'n) = 


e"^+^^Z2L(iV,M) . _ e 


N=0 M=0 

with the normalization factor 


y2L{v,p) 


uN{T])+iiM{-ri) 

T^N,Mi'n) = -^ 7 r(r/) 


(34) 


2L 2L-N 


Z2i(N,M) 


(35) 


N=0 M=0 

are invariant measures for the two-component exclusion process m- 


The second equality in (|34l) follows from the trivial identity e'^'^’''^^l§^^(77) = 

these measures the grandcanonical equilibrium dis¬ 
tributions. The normalization Y 2 L{r',p^), is a homogeneous bivariate Rogers-Szego 
polynomial [ 23 ] and is called the grandcanonical partition function. 

The limits p—?■—ooori /—>■—cxd lead to the pure grandcanonical measures 


qT(^) 


2L 


E 

Ar=0 


e-^C2L(iV) 

X2l{i^) 




2L 


E 

M=0 


e^^C2L{M) 

X 2 M 


7^0,m(^) 


(36) 

(37) 
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with the Rogers-Szego polynomial X 2 l(q!) = From fl^ fol¬ 

lows that 7r^^o(^) ~ 1 s^o(^)^o(^)/^ 2 l(-^, 0) with 7 ^ 0 ( 77 ) = gF;fe=_i,+i( 2 fe-i)afe_ gjnce 
0) = C 2 l{N) one hnds that Q^*{r]) is a prodnct measnre in S^q with 
marginals Q^*{k) = (1 -|- afc(e^g^^“^ — l))/(e‘"g^^“^ -|- 1), reminiscent of the blocking 
measnre of the single-species ASEP on Z [2]. Likewise Q^*{'n) is a prodnct mea¬ 
snre in Sq^ with marginals Q^*{k) = (1 -|- — 1 ))/+ 1). The 

density prohles {ak)^Q and (^fc)o^ in the pnre grandcanonical measnres follow by 
straightforward compntation. One has shock prohles 


( )u,0 

( )o,At 


1 -f 

q 2/c-t-i 

1 -|- e^g“2fe+i 


1 

2 


1 -|- tanh 


- ka ^' 


1 

2 


1 — tanh 


kk- kb\ 

V « J. 


(38) 

(39) 


with the shock width ^ = 1/lng and shock positions Ka = — ^/ln( 3 ')/ 2 , k,b = 

(1 -h u/ \nq)/2. 

In order to describe the self-dnality of the process we dehne for conhgnrations 
77 G the fnnctions 

Q^{r]) = QyiVi) = (40) 


From these fnnctions we constrnct the prodnct 


7V(z) M(z) 


QM ■= n n 

i=l i=l 


(41) 


indexed by z = {x, y}, interpreted as a set of coordinates Xi,yi G A and nnrelated 
to 77 . With this dehnition we are in a position to state the second main resnlt of 
this work. 


Theorem 3.3 Let z and 77 be two configurations of the two-component exclusion 
process defined by m with asymmetry parameter IflR) . The process is self-dual 
with respect to the family of duality functions 

D{z,ri) = Tr-\z)Q^{r]) (42) 

where vr“^(z) is the reversible measure (E^. 

We remark that the reversible measnre (l32l) can be expressed as 

7 r(z) = gZS’[2xi-l-M.,(z)]-E,T^l^^[(2?/.-l-iV,,(z)] ^43^ 

by nsing ffT2l) . Particle nnmber conservation trivially indnces independent dnality 
relations for each combination of particle nnmber pairs (A^, M) = {N{r]), M{r])) and 
{N',M') = (N{z),M(z)) with dnality fnnctions 


( 44 ) 










Therefore we refer to a “family of duality functions” rather than just the “duality 
function”. One has (z, ry) = 0 if iV' > iV or M' > M. Using particle number 

conservation one can construct similar duality functions from := = 

g2Efc=_L+i“fcQ^ and Qy := For the one-component 

ASEP is the duality function of [5] . A family of duality functions for the ASEP 
with second-class particles is given by fH2|) via the replacement 6^ —)■ Ufc. Duality 
and reversibility of the process yield the following corollary, see fll42p : 

Corollary 3.4 For an initial distribution Po{rj) with an arbitrary number N of par¬ 
ticles of type A and M particles of type B we have for the time-dependent expectation 

{Qz{t))p^= ^ (Qz')po^(z;^|z';0) (45) 

^'n,m 

where F{z] t\z'] 0) is the transition probability of the two-component ASEP with N = 
N{z) particles of type A and M = M(z) particles of type B. 

Explicit exact expressions for F(z;t|z';0) have been obtained in [18] for the 
inhnite system. 

Finally we present some simple properties characterizing the invariant measures. 
First we remark that by the dehnition of the process - in which the jumps of the 
A-particles “do not see” whether the neighbouring site is vacant or occupied by a 
fi-particle - expectations of the form ()n m canonical equilibrium 

measure (1^ do not depend on M, i.e., ()nm ~ ■ ■ ■ cik„ )]\fo Likewise, 

{hi-■■ hr, )7v,m = (• • • hr, )o,M does not depend on N. 

The second characterization is a sum rule involving the the canonical invariant 
measures and the duality function. 

Theorem 3.5 Let rjpf be a configuration in with N particles of type A and 
M particles of type B and let z be the coordinate representation of a configuration in 
po,rticles of type A and M' particles of type B. Then for all z G 
and rj G one has the sum rule 

^ ^N,M{'n )Qz{'n ) = Qz'i'n) = ^n,m (^ 6 ) 

with a constant independent of ry and z and canonical stationary distribution 

given by 

We remark that X^'’^' = 0 ii N' > N or M' > M. 
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4 Some tools 


4.1 More notation 

A generic time-dependent probability measure Prob [lit — "H] is denoted by P{r],t) 
or P{'n^). For t = 0 we use the notation Po{rf) := P(r/, 0). If t is irrelevant we omit 
the argument t and write P{'n). We also dehne the transition probability 

:= Prob [r)^ = r/'|r;o = r)] (47) 

from a conhguration r; to a conhguration 7]'. 

The expectation of a function /(r/) is denoted by (/) := J2r) fi'n)P{'n)- If we 
specify time and consider an initial distribution Po(^) we use for the expectation of 
a function f{r]t) the notation 

{fi'nt))po (48) 

V v' 

or simply {fit))p^. For an initial distribution Po{r}') = 5^/,, concentrated on a 
conhguration rj we write ( f{rit) or ( f{t) )^. 

4.2 Matrix form of the generator 

It turns out to be convenient to write the generator (lT5l) in the so-called quantum 
Hamiltonian form [20] which is widely used in the physics literature on stochastic 
interacting particle systems and which was a given a formal probabilistic description 
in [19] . However, this approach does not seem to be well-known in the probabilistic 
literature. For self-containedness and for introduction of our notation we summarize 
the main ingredients. 

4.2.1 Choice of basis, inner product, and tensor product 

In order to write the matrix H explicitly one has to choose an concrete basis, i.e., to 
each conhguration rj one has to assign a specihc canonical basis vector. Following [3] 
we use ternary ordering, i.e., we assign to each conhguration 77 the canonical basis 
vector 

2L 

I'i'n) = ^ + '^v{j - (49) 

i=i 

of the complex vector space C'^ with dimension d = 3^. This basis vector has 
component 1 at position 7 ( 77 ) and 0 else. We work with a vector space over C rather 
than over M since in computations one encounters eigenvectors of H which may be 
complex. 

We denote the basis vectors, which we consider to be column vectors, by 1 77 ). 
We shall also use the notations |z) and |x,y) instead of [ 77 ). The basis vectors 
for conhgurations with a hxed number N of particles of type A and M particles of 
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type B are denoted by | r/^vM )• We define also the dual basis {r]\ = , where 

the superscript T on vectors or matrices denotes transposition. 

The inner product of two arbitrary vectors (w \ with components Wi and (v \ 
with components Vi is defined by 

d 

(tc I n) = ^ WiVi (50) 

i=l 

without complex conjugation. In particular, we have the biorthogonality relation 

( 77 1 77 ') = (5^^, (51) 

Next we introduce the tensor product |n)( 8 )(t(;| = |n)(t(;|. This tensor product 
is a d X d-matrix with matrix elements (| n)(u' \ )ij = ViWj. Specifically we have the 
representation 

1 = (52) 

V 

of the d-dimensional unit matrix, expressing completeness of the basis. 

4.2.2 Measures and expectation values 

A probability measure P{rj) is represented by the probability vector 

|P> = 5^P(77)|7,). (53) 

V 

From the inner product (|5T|) and from (1?I1) we find 

Cf{rj) = -{f\H\'n) (54) 

where the vector ( / | = fi'n){'n I has components /(t 7 ). The semigroup property 
of the Markov chain is reflected in the time-evolution equation 

|Pi) =e-^'|Po) (55) 

of a probability measure Po(t 7 ). 

Normalization implies 

(s|P) = l (56) 

where the summation vector 

(57) 

is the row vector where all components are equal to 1 . As a consequence one has 

(s|P = 0 (58) 

which means that the summation vector is a left eigenvector of H with eigenvector 0. 
This property follows from the fact that a diagonal element of is by construction 
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the sum of all transition rates that appear with negative sign in the same column rj 
of H. The vector corresponding to a stationary distribution is denoted | tt*). This 
is a right eigenvector of H with eigenvalue 0: 

H\n*) = Q. (59) 

and normalization (s | tt* ) = 1. An unnormalized right eigenvector with eigenvalue 
0 is denoted | tt ). 

The expectation (/)p of a function /(r/) with respect to a probability distribu¬ 
tion P{rf) becomes the inner product 

{f)p = {f\P) = {s\f\P) (60) 

where 

f ■■=^f\r))\'n){r)\ (61) 

V 

is a diagonal matrix with diagonal elements /(^y). Notice that 

f{r]) = {rjlflr]) = {s\f\r]). (62) 


For an initial distribution Pq we can now use the dehnitions fl5T]) . fl53l) . fl57l) and the 
representation (l5^ of the unit matrix to recover (l48l) in the matrix form 


(/W)po = 

77 7]' 


71^ 

(63) 

'/ 

= (s|/e-^*|Po) 


Here 

Pi.l'A'n.Q) = (ry'|e“^*|ry) 

(64) 

is the transition probability fl47)). 

For a normalized stationary distribution we also dehne the diagonal matrix 

r := ^7r*(ry)| ry)(ry|. 

(65) 


V 


For ergodic processes with hnite state space one has 0 < vr*( 77 ) < 1 for all ry. Then 
all powers (tt*)" exist. In terms of this diagonal matrix we can write the generator 
of the reversed dynamics as 

Hvev ^ 

Reversibility means = H. An unnormalized stationary distribution vr for which 

Hti = (67) 

holds with 

^ = '^^iv)\'n){r]\. ( 68 ) 

V 

is called a reversible measure. 
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4.2.3 Explicit form of the generator 

In order to write the generator H explicitly we define the following matrices: 


a := 


a := 



&+ ;= 


b- := 



c : = 


c : = 


0 

0 

1 \ 

0 

0 

0 , (69) 

0 

0 

0 / 

0 

0 

0 \ 

0 

0 

0 , (70) 

1 

0 

0 / 


the diagonal projectors 

1 0 0 
a:= \ 0 0 0 
0 0 0 


V := 



b := 



and the three-dimensional unit matrix 

1 = d + i) + b. 


(71) 


(72) 


For matrices M the expression will denote the j-fold tensor product of M 
with itself if j > 1. For j = 1 we define := M and for j = 0 we define = 1 
with the c-number 1. For arbitrary 3 x 3-matrices u we define tensor operators 

Uk := (g) M 0 (73) 

which allow us to write the generator H for the two-component ASEP on the lattice 
{—L -I- 1,..., L} as [3] 


L-l 


H — hk,k+i 

k=—L-\-l 

with the hopping matrices 

dk,k+i ■ ^ (^dkVk-{-i cip, Vkbk-{-i bp. bp_^_^ -|- dkbk-{-i Cp Cp_^_^ 

-\-£ — ttp CLp_^_-p -f- bkVk-\-i — bp bp_^-^ -\- b^dk+i — Cp Cp_^_^'j 

With (lT8l) we split H = Hd + Ho into its offdiagonal part 


(74) 


(75) 


L-l 


Ho = -w [^( 

k=—L-\-l 

and its diagonal part 




Kh+i + Ck c'k+i 


)+^-M 


%+i + h H+i + cuc 


k '^fc+l)] 

(76) 


L-l 


Hd = w 


\q (dkVk+i + Vkbk+1 + dkbk+ij + q ^ y^Vkdk+i + bkVk+i + bkdk+i 


k=—L-\-l 


(77) 

For more details of the construction of H in the tensor basis we refer the reader 

to [3]. 
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4.3 Duality 

We recall the concept of duality in matrix form [211 [10], see also [25] for a detailed 
discussion. In this subsection X and represent arbitrary hnite-dimensional state 
spaces. Consider two processes Xt and ojt and a function D : X xQ C. Notice that 
the function D{x,u) can be understood as a family of functions /j, : hi i—)■ C indexed 
by X and defined by /x(a;) := D{x,u)), or, alternatively as a family of functions 
g^} X C indexed by oj and defined by gu}{x) := D{x,u)). 

The two processes are said to be dual to each other if 

{D{x,ujt))^ = {D{xt,uj))^ (78) 

We remark that with the definitions introduced above we have 

( D{x, ut) )^ = Y1 D{x, 0) = ( f^{t) (79) 

( D{xt, = D{x\ uj)P{x\ t\x, 0) = ( g^{t) (80) 


so that duality can be stated as 


= {gu.it))^ (81) 

with (/x(0))^ = (^..(0))^ = D{x,uj). 

In order to make contact with the quantum Hamiltonian formalism we define 
lx) as a canonical basis vector of and jcj) as a canonical basis vector of 
Let ( s I and ( s | be the corresponding summation vectors. Define the matrix 

D = ^^Il(x,r7)|x)(a;| (82) 

X U) 

with matrix elements {x\D\uj) = D{x,u). The processes Ut and Xt with generators 
H and G are dual to each other w.r.t. the duality function D{x,u) if 

DH = G'^D. (83) 


It is easy to prove the equivalence of this dehnition with the original definition 
(178|) . Since the kind of arguments underlying this equivalence are important for the 
present matrix formulation of duality we present them here in detail: 


{D{x,ujt))^= D{x,u')P{u',t\u,0) (84) 

= (x |Z1| |e“''^*| cu) (85) 

= {x\De~^^\u) ( 86 ) 

= {x\e~^'^^D\uj) (87) 

= (x x')(x'|Z1 | ca) (88) 

= T.X'D{x',uj)P{x',t\x,0) ={D{xt,uj))^ (89) 


14 


In going from (]85|) to (|86|) and from (j88ll to (j8^ we use the representation of the 
unit matrix constructed in analogy to (l52l) . In the step from from fl86|) to flHTD we 
apply the dehnition (18^ . Since we have a chain of equalities, it can be read in both 
directions. Thus the equivalence is established. 

In order to express the alternative dehnition (1811) in matrix form we introduce 


the diagonal matrices 

fx = '^D{x,uj)\uj){uj\, = '^D{x,uj)\x){x\. (90) 

LJ X 

The duality relation fISTll reads 

(s cj) = (s x). (91) 

To prove of equivalence of fIM]) with fl83|) we note that by construction 

{s\fx = {x\D, g^ = \s) = D\uj). (92) 

Then it follows that 

{x\De-^^\u:) (93) 

= {x\e-^^^D\u) (94) 

= (x s) = (s|^^e“^*|a;) (95) 

which establishes the equivalence. 

We end this discussion with a reformulation of Theorem 2.6 of ng. 


Theorem 4.1 Let H be the matrix representation of the generator of an ergodic 
Markov process pt with countable state space and be the matrix form of the 
generator of the reversed process ft- Assume that there exists an intertwiner S such 
that 

SH = (96) 

Then H is self-dual with duality function p) = given by the matrix elements 
of the duality matrix 

D = ^-^S. (97) 

with the diagonal stationary distribution matrix TUR) . 

The proof that SH = implies self-duality with duality matrix D = 

{7r*)~^S is elementary and follows from the chain of equalities 

DH = = H^D. (98) 

The hrst and the third equality are the dehnition fl97)) . the second equality is the 
hypothesis fl9BD of the theorem, and the fourth equality is the reversibility relation 

dnzD. 

Remark 4.2 It follows that if H is reversible then the hypothesis l[9R} reads SH = 
HS, i.e. S is a symmetry of H. Unlike m we do not require S to be invertible. 
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4.4 Representation matrices for u,im] 

4.4.1 Relation between Uq[gl{n)] and Uq[si{n)] 

It is convenient to introduce generators Hj and Hj through 

= L,, R = R - R+i. (99) 

Then the quantum algebra Uq[sl{n)] is the subalgebra generated by 
i = 1,..., n — 1 with relations (|25|) . (l26i) and 

gHb2^-Hi/2 ^ ^-H,/2^H,/2 ^ j 
gHi/2^H,/2 ^ ^Hj/2gHi/2 

(102) 

[X+,X-]=%[H4. (103) 

with the unit I and the Cartan matrix A of simple Lie algebras of type An 

{ 2 * = j 

A,j .= { -1 J = * ± 1 (104) 

[ 0 else. 

The fact that Uq[5\{n)] is a subalgebra of Uq[gl{n)] can be seen by noticing that 
Y^=i R belongs to the center of R[0t„] [21]. 

4.4.2 Tensor representation for n = 3 

In order to distinguish the three-dimensional matrices corresponding to the funda¬ 
mental representation from the abstract generators we use lower case letters. In 
terms of (1^ . (170]) . (ITB the fundamental representation of R[0[(3)] is given by: 

= a^, xf = (105) 

hi = a, h2 = V, hs = b, (106) 

corresponding to 

hi = a — V, h2 = V — b. (107) 

for the representation of the generators Hj of R[s[(3)]. It is convenient to work 
both with hi and the projectors hj expressed in term of the projectors (ITB . 

In terms of the fundamental representation a tensor representation of t/q[s[(3)], 
denoted by boldface capital letters, is given by |3| 

L 

E (108) 

k=—L-\-l 

with 

Y+{k) = q^t-L+Ai-J:tk+Aia+^ (109) 
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and 


with 


Y-{k) 

= Ei=-L+1 

(110) 

y2Hk) 

= qYlt=-L+ik—J2iLk+ikjj— 

(111) 

y2-{k) 


(112) 


L 

H,= ^*(^) 

(113) 


k=—L-\-l 


H^{k) 

= 0 hj (g) 

(114) 


Notice that Hi{k) = and H 2 {k) = Vk — bk- 

For the full quantum algebra 17g[gl(3)] we have the diagonal representation ma¬ 
trices 


Hi= ^k=-.N, H2= ^3= h=-.M. (115) 

k=—L-\-l k=—L-\-l k=—L-\-l 


Here N and M are the particle number operators satisfying 

From these matrices one obtains the representation matrices 

Li = (117) 

The unit I is represented by the 3^-dimensional unit matrix 1 ;= 

The crucial property of the representation 01081) and 0117p that was proved in [5] 
and which is used heavily below are the commutation relations fl25]) which express 
the symmetry of the generator H 07^ under the action of the quantum algebra 
t/j0l(3)]. 


5 Proofs 

5.1 Proof of Theorem (l3.1l) 

(i) We hrst note that uniqueness of follows from ergodicity of the process 

dehned on the subset which is ensured by the fact that the process is a random 
sequence of permutations 

(ii) In [3] we proved, using the quantum algebra symmetry (|28D . that the two- 
component exclusion process dehned by ffT5|) has the unnormalized reversible mea¬ 
sure TT (l32l) . Below we give a direct proof without reference to the quantum algebra 
symmetry. According to the discussion of Section fl4.2p we prove the transformation 
property (1^ with the generator (ITip . Since tt is diagonal one has has Tt~^HdTt = 
for the diagonal part fl77P of H. It remains to show that Tt~^HoTT = Hj for the off- 
diagonal part fl76p . To this end we hrst prove the basic transformation lemma 
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Lemma 5.1 For any finite p ^ 0 we have 

= a± 

^aib 

m L± —a;6 

m _ 

y '^x y y '^x' 


(118) 

(119) 

( 120 ) 
( 121 ) 


Proof: A projector u has the property u = u^. Thus its exponential can be 
written = 1 + (p — l)n. Since aibm is a projector one has 


po^ts+i ^ 1 + (p _ i)ai4tpi. 


( 122 ) 


The tensor construction implies that UkUi = uiUk for fc 7 ^ / and any u. For k = I we 
observe that one obtains by direct computation the relations 


and 


a'^a = b'^a = b~a = c'^a = 0, 
a~v = b~ V = c'^v = c~ V = 0 , 
a'^b = a~b = b'^b = c~b = 0, 

aa~ = ab'^ = ab~ = ac~ = 0, 
-Oa'*" = ab^ = vc^ = vc~ = 0, 
ba'^ = ba~ = bb~ = bc^ = 0, 


a~a = a~, c~a = c~ (123) 

a+b = a'^, b'^v = b'^ (124) 

b~b = b~, c'^b = (125) 

aa+ = a+, ac^ = (126) 

va~ = a~, vb~ = b~ (127) 

bc~ = c~. (128) 


By multilinearity of the tensor product these relations remain valid on each subspace 
k. Relations flllSp - fll2ip then follow from fll22p . □ 

Now we decompose 

TT = Am (129) 

with 


= gSL-L+l(2fc-l)«/=^ = g“EL-L+l(2fc-l)fefe^ fj 


L-1 k 

n n q^ibk-\-i—i>i^k+i (^ 130 ) 

k=—L-\-l 1=1 


Together with one has from (11181) and (11191) of Lemma (I5.ip for 

—L + 1 < k < L 


ia±i-‘= i6±i-‘ = 6f, icfi"* = (131) 

BbfB-' = = aj, Bcjfi-" = (132) 


For the transformation Up one obtains from (11201) and (I12ip of Lemma (I5.ip for 
—L + 1 < k < L 


UafU-^ 


= Q 


T Ei=_L+i i>i^'l2i=k+i 


(133) 
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Ub^U-^ = (-^34) 

Uc^U~^ = q^^tU+iibi-h)±T.tk+i(h-h)(,±_ ('X35) 

Putting these results together and using the projector property fll22p together 
with fll23p - fll28p yield for —L + l<k<L — 1 

(136) 

(137) 

and 


Uatat+iU ‘ = ‘ = C7cfcf+,C/ ‘ = 


(138) 


Since (a^)^ = (and similarly for 6^ and c^) applying the decomposition fll29p 
to the individual terms in flTip yields = Hj and therefore reversibility of tt. 

(iii); We complete the proof of Theorem fl3.ip by proving the normalization 
factor. For a conhguration z = {x, y} dehne yi == yi — Ny.{z). We have by 
dehnition of the partition function 


Z2l{N,M)= ^7r(zM,M)= 


E£i(2xi-i)-E"i(2yi+^-i) 


(139) 


^N,M 


^N,M 


Consider the points r in the Weyl alcove hF|h = —L < xi < ... < xk < T}. We 
also dehne the punctuated Weyl alcove hF|'^(r) = \r for r G hF|h. This allows 
us to write = E.-ew- ^yewl^H^Y 

Next observe that by construction Y.y&wjj^ix) fiVi) = /(l/i)- Therefore 

Z2l{N,M)= Y (140) 


which implies that Z 2 l{N,M) = Z 2 l{N, 0 )Z 2 l-n{ 0 ,M). a classical result from 
the theory of integer partitions [26] yields for the single-species partition functions 
Z 2 l{N, 0) = C 2 L(iV), Z2 l(0, M)C 2 l(M) with the g-binomial coefhcient Ck{N). Ob¬ 
serving that C 2 l{N)C 2 l-n{M) = C 2 L{Af,M) concludes the proof. □ 


5.2 Proof of Theorem (l3.3l) 

5.2.1 Reformulation of the problem 

Step 1: We hrst apply the general considerations of Sec. fl4.3l) to the present case 
of the two-component ASEP. It is convenient to use the occupation variable pre¬ 
sentation Tf^ for one process and the coordinate representation z^ for the dual. The 
duality function, given by (z |Z1| 77 ) in terms of the duality matrix D, is therefore 
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denoted by D{z,r}). If self-duality is valid for some duality function D{z,r}) then 
according to fl90l) we can define a diagonal matrix such that 

D{z,r]) = {s\D^\r]) (141) 

with the summation vector (s |. Then self-duality yields 

{s\D^e~^^\r]) = ^(z'|e“^*|z)(s|T)z|77) (142) 

and, as a consequence from reversibility, Corollary fl45|) . 

Step 2: In [3] we have established the symmetry of the generator under the 
action of f/q[g[(3)]. Moreover, we have reversibility H = of the two-component 
ASEP with the reversible measure 032 p . Then for any matrix S satisfying [S, H] = 0 
Theorem 04. ip yields a duality function 

D{z,r]) = (z|(7 r*(z))"^^|77) = 7 r"^(z)(z |^| r/). (143) 


which means that we can construct duality functions from the symmetry operators 
of the model, i.e., from the tensor representation OlOSp . 0113p . 

Step 3: On the other hand, from 014ip . one has D{z,r]) = {s\Dz\r]) for all 
r] G Therefore one can express the duality function 

D{z,r]) = 7i-^{z)Q^{r]) = 7i~^{z){s\Q^\'n). (144) 


of Theorem 03.3p in terms of a diagonal matrix Qz satisfying 


{z\S = {s\Q^ 


(145) 


and {s\Qz\rj) = Qz{r}) given in OTTl) . Therefore the task at hand is to hnd a 
symmetry operator S that satishes (I145p with the diagonal matrix Qz with matrix 
elements given by (ITTll . 

Step 4: In order to choose S we observe that D{(/),ri) = 1, corresponding to 
Dq = = 1. The non-trivial information one gains is that (0|*S' = (s| which 

means that the symmetry operator S generates the summation vector from the 
vacuum vector (0 |. From the explicit representation obtained in [3] we hnd as a 
candidate 


2 ^ 2^ MU ImlJ 


(146) 


Since 7 r(z) is known the remaining task is to construct D{z,r]) as stated in the 
theorem by proving (I145p . 


5.2.2 Technical lemmas 

We prove the following lemmas: 
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Lemma 5.2 Consider coordinate sets x' = x U r and y' = y U s. For k ^ r and 
I ^ s one has 


A^fc({xUr, •}) = A^fcdx, •}) + Arfc({r, •}) (147) 

= ]Vi({x,.}) + ]V({r,.})- Yi eftn) (148) 

i=l 

Mi{{;yUs}) = Mi{{;y}) + Mi{{;s}) (149) 

M({.,s}) 

= MK{-,y}) + M({-,s})- Q{l,Si). (150) 

i=l 

Proof: The function Q{r,x) defined in ([6]) satisfies 

0(r, x) = 1 — 0(x, r) — 6r,x (151) 

x—1 L 

dr,k = C){r,x), 5r,k = C){x,r) (152) 

k=—L-\-l k=x-\-l 

From (fTOl) . (|T^ and (fT3|) we have 

V(z) L N{z) 

Y = ^(^) “ ^'■(^) “ Y 

i=l l=r+l i=l 

M(z) L M(z) 

Y Y1 ~ ~ 

i=l l=r+l i=l 

Specihcally, for iV(z) = 1 or M(z) = 1 resp. we obtain from 01511) - 01541) 

l\^r-({2:, •}) = Mr{y) = Q{y,r), (155) 

and more generally one hnds for z = {x, y} from 0131) and 0152p 

7V(z) M(z) 

A^^({x, •}) = ^0(xi,r), Mr{{-,y}) = Y^iViC)- (156) 

i=l i=l 

With 01551) and 0156p one then hnds 

N{z) M(z) 

W({x, •}) = W({a:i,-}), Mr{{-,y}) = Y^ri{-,yi})- (157) 

i=l i=l 

The hrst equality 01471) in the lemma then follows from 01571) . The second 
equality 01481) arises from 015ip and 01561) . bearing in mind that by assumption 
/c ^ r. The proof of 0149p and 0150p is analogous. □ 
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In particular, for x = y = 
the inversion formulas 


■}) 




0 one obtains from Lemma (I5.2p for k ^ r and I ^ s 


ivHi-.-})- E 

2=1 

(158) 


(159) 


2=1 


We also derive the projector lemma. 

Lemma 5.3 The tensor occupation operators hk, act as projectors 

N{r,) 

41»?) = Ofcl r/) = ^ 77 ) (160) 

2=1 
M{ri) 

hlr]) = bk\'n) = '^ 6 y^,k\'n) ( 161 ) 

2=1 

with the occupation variables and b^ m (or particle coordinates x* and Ui respec¬ 
tively) understood as functions of rj or z = 77 . 

Proof: The first equality in each equation is inherited from the definition of the 
projectors fITTll by multilinearity of the tensor product, the second equality follows 
from flT^ . □ 

Finally we note two combinatorial identities for sums over the permutation group 
Sn- One has 

-(0(r-„rP)+n(n-l)/2 _ ^-2^"=! Et? 

cr€Sn 

which can be proved by induction using [n]g = X]fc=o 

L L 

E E /(ri,...,r„) = EE /(a(ri,...,r„)) (163) 

ri=—L+l rn=—L+l r„ 

for functions that vanish whenever r^ = Vj. Here the sum over r„ denotes the 
summation over the Weyl alcove 


5.2.3 Main steps 

After these preparations we go on to prove for z = {x, y} the property 

W({x,y}) M( x,y}) 


!x,y4={si n n ^ 


(164) 


2=1 


2 = 1 
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The matrices on the r.h.s. are the operator form of the functions Q]^ ^ dehned in 
fHOj) . According to fll45p proving (11641) proves the theorem. 

From the representation (IllOl) . from (11601) and from the dehnition (IT4|) one hnds 

(x, y |Fi“(r) = U r, y |. (165) 

By iteration 

(x,y|y'f(ri)...yi"(r„) = x U r„, y | (166) 

with the dehnition fq = 0. Since from ffT4)) . fll47p . fll5ip one can write 

i-i 

Ar. ({x U rj_i, •}) = 2Nr. ({x, •}) + 2 ^ ©(r^, rj) - (iV(x, •}) + j - 1) (167) 

i=l 


one has 


n j-1 


^A^,({xUrj_i,-}) = 2^iV^^({x,-})-niV({x,-})+2^^0(ri,rj)--n(n-l). 


i=i 


j=i 


j=i i=i 


(168) 

Now we observe that the term ({x, •}) — iV({x, •})) is invariant under 

permutations of the coordinates rj. Using the fact that {Y^{r)y = 0 and the 
combinatorial properties (11621) and (11631) then yields 


< -. y 1^ = E 9" (3, u r„,: 

mJ'?' 


(169) 


The next step is to invoke Lemma (15.21) to express ({x, •}) in terms of single¬ 
particle step functions Nx^{{rj, ■}) with inverted arguments. This initiates the fol¬ 
lowing chain of equalities for the exponent F'r({x, •}) := — ({x, •}) — 

^({x, •})] of q: 

n / ^({x,-}) 

Sr({x,.}) = -5 ^(jv({x..})- 2 ^ lV,,({r^,.}) 

i=l 


n Af({x,-}) 

(21V,,({r„.})-!) 

j=l i=l 

i=l 


Since trivially 


{xUr„,y I ^ {xUr„,y .. .a 


(170) 


xn 


(171) 
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one arrives at 


(x,y 



'Af({x,y}) 

I n I (xUr„,y 

i=l 


'Af({x,y}) 


^^xi({r,y}) 


^ n I Un,y 


\ *=1 


( 172 ) 


where the summation in the second equality has been changed to the extended 
Weyl alcove with n = iV({x, y}) +n. This is possible since the product of indicators 
dxi ■ ■ ■ Oxjv cancels all terms not belonging to the original Weyl alcove n. 

Next one uses the dehnitions m, ma and the projector property (11601) to 
express y43;.({r,y}) in terms of projectors. The result is 


,x,y 


(n- 


W- 


/^({x,y}) 

5^(rn,y| n ^ 

bi \ *=i 

Af({x,y}) ^ 

^(rn,y| n 

ff,. 7=1 


E Xj 1 -s V ^ 

i^l ^j-z2j=Xi + 1^3 ^ 


Xi 


(173) 


Using the commutator relation [T^ , ] = 0 (IMj) and going through similar 

steps yields 


,x,y 


(Uf)- (y+y 


[n\g\ [m\g\ 


W({x,y}) M({x,y}) 


y~! y~! ( ^ri, Sm I j Qxi Q 


(174) 


' n 


2=1 


2=1 


with m = M ({x, y}) + m and 


Qy = q (175) 

Observe now that the summation on the r.h.s involves only the vector (r^i, |. 
Since the summation is over the Weyl alcove one has 

EE ( ^hi I ( '^r2,m |* (^^6) 

^ n ^rh 


Therefore 


(x,y 


(y-)n 
[n]g\ [m]g\ 


^n.m. 


'7V({x,y}) M( x,y}) 

If n Qn 

2=1 2=1 


(177) 


Notice that on the r.h.s. the only dependence on the n and m is in the summation 
vector {sn^m I for the sector with n = iV({x, y}) + n particles of type A and fh = 
M({x, y}) + m particles of type B. 


24 









The final step is to take the double sum (11461) . Terms such that n + m > 
2L - iV({x,y}) - M({x,y}) are zero since {sN, 2 L-N\yr = {sN, 2 L-N\y'^ = 0, 
corresponding to the exclusion principle that forbids creating configurations with 
more than 2L particles on A. On the other hand, 

/7V({x,y}) M({x,y}) \ 

(On.™ I n Qi. n Qy.] = 0 if n < A^({x, y}) or m < M({x, y}), (178) 
\ i=l i=l I 

due to the projectors contained in the operators This yields fll64p with 

Qz = QxQy (179) 

and proves Theorem fl3.3p by taking the scalar product with { s \ and 1 77 ). □ 

5.3 Proof of Theorem (13.51) 

The first equality follows from duality and ergodicity by taking the limit t —?■ 00 
in the expectation fl45l) with coordinate sets z' representing configurations with N' 
particles of type A and M' particles of type B. 

We also have from fl45l) by taking Pq = ^y considering the canonical 

equilibrium distribution for N particles of type A and M particles of type B, 


= (z|e-^*|V)v,M). ( 180 ) 

^N,M v I I i / 

with a vector 

\Vn,m) = ( 181 ) 

that has support in the subspace corresponding to configurations with N' particles 
of type A and M' particles of type B. 

Because of stationarity the l.h.s. does not depend on time. This implies by 
ergodicity that | Vat^m ) on the r.h.s. must be proportional to the (unique) stationary 
probability vector for configurations in Hence 

~ ^N,M ^ (182) 

where is some constant with N' = N{z) = N{z') and M' = M(z) = M(z'). 

This proves the second equality in the theorem. □ 


Acknowledgements 

This work was supported by DFG and by CNPq through the grant 307347/2013-3. 
GMS thanks the University of Sao Paulo, where part of this work was done, for kind 
hospitality. 


25 











References 


[1] Ferrari, P.A., Kipnis, C., Saada, E.: Microscopic Structure of Travelling Waves 
in the Asymmetric Simple Exclusion Process. Ann. Probab., 19(1) 226-244 
(1991) 

[2] Liggett, T.M.: Interacting particle systems. Springer, Berlin, (1985) 

[3] V. Belitsky and G.M. Schiitz, Quantum algebra symmetry and reversible mea¬ 
sures for the ASEP with second-class particles, (2015) 

[4] Schiitz, G., Sandow, S.: Non-abelian symmetries of stochastic processes: deriva¬ 
tion of correlation functions for random vertex models and disordered interact¬ 
ing many-particle systems. Phys. Rev. E 49, 2726-2744 (1994) 

[5] Schiitz, G.M.: Duality relations for asymmetric exclusion processes. J. Stat. 
Phys. 86(5/6), 1265-1287 (1997) 

[6] Belitsky, V., Schiitz, G.M.: Diffusion and coalescence of shocks in the partially 
asymmetric exclusion process. Electron. J. Prob. 7, Paper No. 11, 1-21 (2002) 

[7] Imamura, T., Sasamoto,T.: Gurrent Moments of ID ASEP by Duality. J. Stat. 
Phys. 142(5), 919-930 (2011) 

[8] P. Le Doussal and P. Galabrese, The KPZ equation with flat initial condition 
and the directed polymer with one free end, J. Stat. Mech. (2012) P06001. 

[9] Giardina, G., Kurchan, J., Redig, F.: Duality and exact correlations for a model 
of heat conduction. J. Math. Phys. 48(3), 033301 (2007) 

[10] Giardina, G., Kurchan, J., Redig, F., Vafayi, K.: Duality and Hidden Symme¬ 
tries in Interacting Particle Systems. J. Stat. Phys. 135, 25-55 (2009) 

[11] J. Ohkubo, Duality in interacting particle systems and boson representation, 
J. Stat. Phys. 139 454-465 (2010). 

[12] Garinci, G., Giardina, G., Giberti, G., Redig, F.: Duality for Stochastic Models 
of Transport J. Stat. Phys. 152(4), 657-697 (2013) 

[13] Garinci, G., Giardina, G., Giberti, G., Redig, F.; Dualities in population ge¬ 
netics: A fresh look with new dualities. Stochastic Processes Appl. 125(3), 
941-969 (2015). 

[14] G. Garinci, G. Giardina, F. Redig, and T. Sasamoto. A generalized Asymmet¬ 
ric Exclusion Process with Uq(sl2 ) stochastic duality. arXiv: 1407.3367 2014. 
http://arxiv.org/abs/1407.3367, Gited 13 Mar 2015. 

[15] Borodin, A., Gorwin, L, Sasamoto, T.: From Duality to Determinants for Q- 
TASEP and ASEP. Ann. Probab. 42, 2314-2382 (2014). 


26 


[16] I. Corwin and L. Petrov, Stochastic higher spin vertex models on the line. 
arXiv:1502.07374 [math.PR] (2015) http://arxiv.org/abs/1502.07374, Cited 13 
Mar 2015. 

[17] S. Chatterjee, G. M. Schiitz Determinant representation for some transition 
probabilities in the TASEP with second class particles, J. Stat. Phys. 140, 
900-916 (2010). 

[18] C. A. Tracy, H. Widom, On the asymmetric simple exclusion process with 
multiple species, J. Stat. Phys. 150, 457-470 (2013). 

[19] Lloyd P, Sudbury A, and Donnelly P 1996, Quantum operators in classical 
probability theory: I. “Quantum spin” techniques and the exclusion model of 
diffusion. Stock. Processes Appl. 61(2), 205-221. 

[20] G. M. Schiitz, Exactly solvable models for many-body systems far from equilib¬ 
rium, Phase Transitions and Critical Phenomena. Vol. 19, Eds. C. Domb and 
J. Lebowitz, Academic Press, London, (2001). 

[21] M. Jimbo, A q-Analogue of U{gl{N -|- 1)), Hecke Algebra, and the Yang-Baxter 
Equation, Lett. Math. Phys. 11, 247-252, 1986 

[22] C. Burdik, M. Havhcek, and A. Vancura, Irreducible Highest Weight Repre¬ 
sentations of Quantum Groups U{gl{n,C)) Comm. Math. Phys. 148(2) (1992), 
417-423. 

[23] Gasper, G., Rahman, M.: Basic hypergeometric series. Encyclopedia of Math, 
and its Appl. 96 (2nd ed.), Gambridge, Gambridge University Press (2004) 

[24] A. Sudbury and P. Lloyd. Quantum operators in classical probability theory. 
11: The concept of duality in interacting particle systems. Ann. Probab. 23(4) 
(1995), 1816-1830. 

[25] S. Jansen and N. Kurt, On the notion(s) of duality for Markov processes. Prob. 
Surveys 11 (2014), 59-120. 

[26] Andrews, G.E.: The Theory of Partitions, Encyclopedia of Math, and its Appl. 
2, London, Addison Wesley (1976) 


27 


